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Subject:  Correcting  for  the  Sidewall  Boundary  Layer  in  Subsonic 
Two-Dimensional  Airf oil/Hydrofoil  Testing** 

Abstract:  Historically,  two-dimensional  airfoil  or  hydrofoil  section 

characteristics  have  been  obtained  by  measuring  individually 
the  lift,  drag  and  pitching  moment  by  the  most  accurate 
technique  available.  The  use  of  force  balances  to  measure 
the  three  quantities  simultaneously  has  met  with  only 
partial  success.  Although  the  lift  and  pitching  moment 
data  have  usually  been  acceptable,  the  drag  data  have 
varied  by  as  much  as  an  order  of  magnitude  from  the  accepted 
NACA  reference  data.  To  investigate  the  parameters  which 
influence  two-dimensional  force  measurement  and  force  balance 
design,  an  experimental  program  was  conducted  in  the  subsonic 
wind  tunnel  of  the  Applied  Research  Laboratory  at  the 
Pennsylvania  State  University  (ARL/PSU) .  From  the  results 
of  this  test  program  the  sidewall  boundary  layer  was  identified 
as  the  primary  factor  contributing  to  the  erroneous  drag 
measurements.  A  correction  procedure  which  is  based  on  the 
airfoil/hydrofoil  geometry,  the  flow  environment  and  the 
measured  data  was  developed.  Corrected  data  from  the  subject 
test  program  and  from  similar  programs  in  other  experimental 
facilities  for  both  symmetrical  and  cambered  sections  are  in 
good  agreement  with  NACA  data  in  all  cases. 


*Aerodynamicist,  Goodyear  Aerospace  Corp.,  Akron,  OH,  formerly 
Graduate  Assistant,  ARL/PSU 


**This  paper  is  being  prepared  for  presentation  at  the  AIAA  12th 
Aerodynamic  Testing  Conference  in  Williamsburg,  VA,  March  1982. 
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NOMENCLATURE 

2 

aspect  ratio  =  s  /sc 
exponent  in  Equation  (11) 
airfoil  chord  length 

sectional  drag  coefficient  =  D/Cq^  •  s) 

sectional  drag  coefficient  at  c^  =  0.0 

sectional  lift  coefficient  =  L/(qM  •  s) 

local  slope  of  the  c^  versus  a  curve 

slope  of  the  linear  portion  of  the  c^  versus  a  curve 
(usually  evaluated  in  the  c^  =  0.0  region) 

the  drag  force 

Hawthorne's  approximation  of  the  energy  in  a  secondary  flow 
(Equation  (1)) 

defined  by  Equation  (2)  and  Figure  12 

functional  operators  used  in  developing  Equation  (12) 

proportionality  constants  used  in  developing  Equation  (12) 

the  lift  force 

defined  by  Equation  (3) 

atmospheric  pressure 

static  pressure 

total  pressure 

2 

local  dynamic  pressure  =  l/2pu 

2 

free  stream  dynamic  pressure  =  l/2pVro 
Reynolds  number  =  V^c/v 
airfoil  span 

maximum  airfoil  thickness 

local  velocity 

free  stream  velocity 
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NOMENCLATURE  (Cont'd) 

x  =  distance  parallel  to  test  section  centerline  measured  from 

the  leading  edge  of  the  two-dimensional  test  chamber 

a  =  airfoil  angle  of  attack 

ACj  =  required  correction  in  c^ 

6  =  boundary  layer  thickness  at  the  balance  shaft  location 

* 

6  =  boundary  layer  displacement  thickness  at  the  balance  shaft  location 

p  =  mass  density  of  the  fluid 

v  =  kinematic  viscosity  of  the  fluid 
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INTRODUCTION 


In  the  early  1970’s  the  Applied  Research  Laboratory  at  the  Pennsylvania 
State  University,  in  an  effort  to  develop  new  propulsor  blade  design  criteria, 
attempted  to  measure  the  lift,  drag  and  pitching  moment  characteristics  of  a 
hydrofoil  that  spanned  the  two-dimensional  test  section  of  the  ARL/PSU  12-inch 
(304.8  mm)  cavitation  tunnel.  The  measurements  were  conducted  with  a  three- 
component,  cantilevered,  mechanical  force  balance*.  With  the  exception  of  the 
pitching  moment  characteristics,  the  data  were  in  disagreement  with  accepted 
National  Advisory  Committee  for  Aeronautics  (NACA)  section  data  and  with 
measurements  published  by  other  investigators.  In  an  effort  to  resolve  these 
differences  and  to  develop  a  means  of  accurately  measuring  the  desired 
hydrofoil  force  and  moment  parameters,  a  basic  research  program  which  was 
conducted  by  Jacobs  [1]**  in  the  ARL/PSU  subsonic  wind  tunnel  was  initiated. 
This  facility  provided  an  environment  in  which  the  aerodynamic  and  geometric 
parameters  common  to  both  airfoil  and  hydrofoil  testing  could  be  easily 
and  economically  varied. 

Although  air  and  water  are  both  fluids,  fundamental  differences  exist 
in  the  testing  requirements  between  the  two  media.  Testing  airfoil  shapes 
in  water  (hydrofoils)  introduces  additional  problems  such  as  the  handling 
of  larger  gross  forces  and  waterproofing  requirements.  However,  the  major 
concern  in  water  tunnel  measurements  is  the  cavitation  phenomenon.  A  hydro¬ 
foil  may  operate  in  any  or  all  of  three  different  flow  regimes;  namely, 
fully  wetted  flow,  partially  cavitating  flow  or  fully  cavitating  flow. 

Testing  in  fully  wetted  flow  differs  little  from  low-speed  wind  tunnel  testing, 
and  it  is  this  regime  that  was  of  primary  concern  in  the  original  ARL/PSU 
water  tunnel  test  program. 

A  review  of  the  literature  indicated  that  most  of  the  currently  used 
NACA  airfoil  section  data  were  measured  in  air  tunnels  during  the  1930' s 
and  1940's  at  both  the  Langley  and  Ames  Research  Centers.  The  majority  of 
these  data  were  acquired  by  measuring  individually  the  lift,  drag  and 
pitching  moment  by  the  most  accurate  means  available.  Generally,  this 
meant  either  the  measurement  of  the  pressure  distributions  on  the  ceiling 
and  floor  of  the  test  section  or  the  use  of  a  force  balance  to  obtain 
lift.  Drag  data  were  obtained  from  either  wake  surveys  or  via  surface 
pressure  distributions;  a  torsional  balance  was  usually  used  to  measure 
the  pitching  moment.  Typical  of  these  measurement  programs  were  those 
conducted  by  Loftin  and  Smith  [2]  in  1949. 

In  a  water  tunnel,  placing  pressure  taps  on  a  hydrofoil  surface 
or  on  the  test  section  walls  can  cause  premature  cavitation  and  result 
in  erroneous  pressure  measurements.  At  certain  flow  conditions,  pressure 
probes  in  the  profile  wake  are  also  subject  to  cavitation  problems. 

For  these  reasons  water  tunnel  force  measurements  are  best  performed 
by  a  mechanical  force  balance.  With  such  a  balance,  forces  can  be 
measured  directly  without  marring  the  models  surface.  Mechanical  balances 
are,  however,  not  totally  free  of  problems.  Balances  measure  all  forces 
applied  to  a  model.  If  forces  occur  on  a  model  which  are  not  those 
associated  with  two-dimensional  flow  (due  to  such  causes  as  secondary 


*In  this  paper  a  cantilevered  balance  refers  to  the  experimental  configura¬ 
tion  where  the  model  is  supported  from  only  one  side  of  the  test  section 
and  the  supporting  shaft  is  equipped  with  a  force  measuring  device. 

**Numbers  in  brackets  refer  to  references  at  end  of  report. 
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flows  where  the  model  intersects  a  tunnel  wall)  the  balance  will  also 
measure  them.  A  correction  procedure  is  then  required  to  reduce  this 
measured  total  force  to  a  two-dimensional  force.  Despite  the  required 
corrections,  a  mechanical  balance  appears  to  be  the  best  method  for 
measuring  forces  on  hydrofoils  and  was  selected  as  the  primary  measuring 
instrument  for  the  ARL/PSU  wind  tunnel  investigation. 

Other  investigators  have  used  mechanical  balances  for  hydrofoil  force 
measurements.  Kermeen  [3]  used  a  cantilevered  balance  system  at  the 
California  Institute  of  Technology  (CIT)  to  measure  forces  on  NACA  4412 
and  Walchner  profile  7  hydrofoils.  In  correcting  his  data,  Kermeen 
included  the  effect  of  tare  forces  and  four  of  the  traditional  sources 
of  error  described  by  Allen  and  Vincenti  [4];  namely,  solid  blockage, 
wake  blockage,  lift  effect  and  horizontal  buoyancy.  Kermeen' s  measured 
lift  and  pitching  moment  data  for  the  NACA  4412  hydrofoil  were  in  good 
agreement  with  previous  NACA  data.  The  drag  data  agreed  well  in  the 
lower  drag  range,  but  for  angles  of  attack  greater  than  five  degrees 
Kermeen' s  results  were  nearly  a  factor  of  two  higher  than  the  NACA 
data  as  shown  in  Figure  1. 

Earlier  measurements  of  a  similar  nature  were  conducted  by  Daily  [5] 
at  CIT.  Daily  used  a  cantilevered,  three-component  balance  which  was  a 
forerunner  of  Kermeen' s  balance  to  measure  the  forces  on  a  NACA  4412 
hydrofoil  in  the  High-Speed  Water  Tunnel  (HSWT) .  A  full-span  hydrofoil 
was  used  at  low  angles  of  attack;  to  halve  the  forces  at  higher  angles 
of  attack  a  split  foil  was  used.  Daily  implied  that  his  data  required 
little  correction.  As  with  Kermeen' s  data  the  lift  and  pitching  moment 
measurements  were  in  good  agreement  with  the  NACA  data.  Daily  was 
primarily  interested  in  the  low  angle  of  attack  range  where  the  drag 
measurements  were  in  good  agreement  with  the  NACA  data.  However,  at 
higher  angles  of  attack  his  drag  measurements  were  approximately  50.0 
percent  high.  Figure  2. 

Since  measurement  of  the  pitching  moment  had  not  been  a  problem 
in  the  past,  a  two-component  cantilevered  balance  was  designed  and 
fabricated  for  the  ARL/PSU  wind  tunnel  test  program.  The  resulting 
measurements  [1]  produced  lift  data  (with  the  traditional  corrections 
applied)  that  were  in  agreement  with  NACA  data.  Figure  3.  However, 
the  corrected  drag  data.  Figure  4,  show  differences  similar  to  those 
reported  by  Kermeen  and  Daily.  Subsequent  pressure  surveys  at  the 
midspan  of  the  airfoil  resulted  in  drag  data  that  were  in  agreement 
with  the  reference  NACA  two-dimensional  section  data.  Figure  5.  After 
experimental  study  of  several  error  sources,  it  was  reasoned  that  the 
two-dimensionality  of  the  flow  was  being  altered  by  the  presence  of  the 
sidewall  boundary  layer. 

In  wind  tunnel  testing  the  sidewall  boundary  layer  can  be  controlled 
by  blowing  or  suction  techniques  to  yield  the  desired  two-dimensional 
data.  Because  these  techniques  can  aggravate  the  cavitation  problem  In 
hydrofoil  testing,  the  alternate  approach  of  using  a  correction  procedure 
is  preferred. 


-7- 


25  August  1981 
ALT : PPJ : GBG : cag 


The  independent  work  of  Barber  [6]  and  Hawthorne  [7]  relative  to 
the  secondary  flows  at  a  strut-wall  intersection  as  a  function  of  the 
incoming  wall  boundary  layer  characteristics  was  a  starting  point  for 
the  development  of  the  desired  correction  procedure.  Their  work,  when 
combined  with  an  empirical  parametric  evaluation  of  the  Kermeen,  Daily, 
and  ARL/PSU  balance-measured  drag  data  resulted  in  a  correction  procedure 
for  the  effects  of  the  sidewall  boundary  layer  that  was  developed  by 
Jacobs  [1]  and  is  summarized  here  in  Equation  (12).  When  this  correction 
procedure  is  applied  to  the  three  previously  discussed  sets  of  drag  data 
and  to  additional  drag  data  measured  by  Ward  [8]  on  a  NACA  16-309  hydrofoil 
at  C1T,  the  results  are  in  good  agreement  with  NACA  two-dimensional  section 
characteristics,  Figures  13,  14,  and  17.  Thus,  based  on  the  rather  limited 
comparative  data,  it  appears  that  a  correction  procedure  to  account  for 
the  effects  of  the  sidewall  boundary  layer  on  balanced-measured  airfoil/hydro¬ 
foil  drag  data  has  been  formulated. 

The  specific  details  of  the  entire  ARL/PSU  wind  tunnel  test  program 
are  documented  in  the  report  by  Jabobs  [1].  Discussed  in  the  remainder 
of  this  paper  are  the  portions  of  the  ARL/PSU  test  program  relevant  to 
the  lift  and  drag  measurements  by  the  cantilevered  balance,  the  develop¬ 
ment  of  the  correction  procedure  and  the  application  of  the  correction 
procedure  to  existing  experimental  data. 

ARL/PSU  WIND  TUNNEL  TEST  PROGRAM 


Test  Facility  and  Experimental  Hardware 

The  ARL/PSU  subsonic  wind  tunnel  in  which  this  test  program  was 
conducted  is  shown  schematically  in  Figure  6  and  described  in  detail 
in  Reference  [9].  This  facility  is  a  closed  circuit,  closed  jet  air 
tunnel  with  a  octagonal  test  section  which  is  4.0  ft  (1.219  m)  across 
the  flats  and  is  16.0  ft  (4.877  ra)  long.  The  test  section  velocity  can 
be  varied  continuously  up  to  120.0  fps  (36.576  m/sec).  Honeycomb  and 
screens  used  in  the  settling  section  reduce  the  turbulence  level  in 
the  test  section  to  0.06  percent  of  the  free-stream  velocity  at  80.0  fps 
(24.384  m/sec).  For  this  test  program  two  4.0  ft  x  8.0  ft  (1.219  m  x  2.438  m) 
wooden  panels  were  mounted  vertically  18.375  in  (466.725  mm)  apart  to  create 
a  two-dimensional  test  section  as  shown  in  the  balance  installation  drawing. 
Figure  7. 

Two  airfoil  models  of  different  aspect  ratio,  AR,  were  fashioned 
for  use  with  the  cantilever  balance.  A  NACA  0012  airfoil  was  machined 
from  aluminum.  This  model  had  a  chord  of  9.0  in  (228.6  mm)  and  span 
of  18.375  in  (466.7  mm)  for  an  AR  of  2.04.  A  second  model  was  made  of  a 
two  part  expandable  urethane  foam  (span  =  18.375  in  (466.7  mm); 
chord  =  18  in  (457.2  mm);  AR  =  1.02).  When  installed  in  the  test  chamber 
the  airfoil  was  attached  at  its  midchord  to  the  cantilever  balance.  The 
balance  shaft  was  located  midway  between  the  floor  and  ceiling  of  the  test 
section  and  28.0  in  (711.2  mm)  downstream  from  the  leading  edge  of  the 
wooden  panels.  Provisions  were  made  to  mount  the  9.0  in  (228.6)  airfoil 
directly  to  the  balance  shaft  or  mounted  to  an  11.0  in  (279.  4  mm)  diameter 
disk  on  the  balance  shaft. 


-8- 


25  August  1981 
ALT : PP J : GBG : cag 


Ins  truraentation 


The  lift  and  drag  forces  were  measured  by  a  two-component,  canti¬ 
levered  balance  that  is  sketched  in  Figure  7.  This  balance  used  the 
compact  strain-gaged  tension-member  concept  developed  by  Gurney  [10]. 

The  balance  rotated  with  the  airfoil  and  sensed  forces  normal  to  and 
parallel  with  the  chordline.  The  corresponding  lift  and  drag  values 
were  computed  via  standard  vector  resolution  equations. 

For  the  force  measurements  the  reference  velocity,  V  ,  was  measured 
by  a  0.25  in  (6.35  mm)  diameter  pitot-static  probe  that  was  located  on 
the  tunnel  floor  midway  between  the  sidewalls  and  in  the  same  vertical 
plane  as  the  midchord  of  the  airfoil.  The  probe  tip  was  above  the  floor 
boundary  layer. 

Wake  traverses  to  evaluate  the  drag  coefficient  were  conducted  with 
a  0.125  in  (3.175  mm)  diameter  kiel  probe  that  was  located  midway  between 
the  sidewalls  and  in  a  plane  one  chordlengf.h  downstream  from  the  airfoils 
trailing  edge.  For  these  tests  the  reference  pitot-static  probe  was 
located  in  the  traverse  plane  midway  between  the  floor  and  the  centerline 
of  the  test  section.  The  same  kiel  probe  was  used  to  make  sidewall  boundary 
layer  measurements  at  the  balance  shaft  location.  For  these  measurements 
the  reference  pitot-static  probe  in  the  tunnel  floor  was  used. 

To  obtain  horizontal  buoyancy  corrections,  the  sidewall  static  pressure 
gradient  was  measured  by  four  .032  in  (.813  mm)  diameter  static  pressure 
taps.  Beginning  1.0  ft  (304.8  mm)  downstream  of  the  sidewalls  leading 
edge  and  spaced  at  2.0  ft  (609.6  mm)  intervals  thereafter.  The  static 
pressure  taps  were  located  1.0  ft  (304.8  mm)  above  the  test  section  floor. 

The  angle  of  attack,  a,  was  measured  with  a  gunner's  quadrant  which 
was  used  in  conjunction  with  an  airfoil  template.  The  template  was 
machined  so  that  when  the  template  made  a  three-point  contact  with  the 
airfoil's  upper  surface  the  top  of  the  template  was  parallel  to  the 
chordline.  This  technique  enabled  the  angle  of  attack  to  be  measured 
within  ±0.056° . 

Based  on  a  Student-t  analysis  the  greatest  error  in  the  balance- 
measure  c^  values  was  ±0.0038  and  in  c^  was  ±0.0007.  A  similar  analysis 
for  the  c^  values  obtained  by  the  wake  traverses  indicated  a  greatest 
error  of  ±0.00014. 

Measurements  and  Results 


The  flow  characteristics  of  the  two-dimensional  test  section  were 
established  by  a  series  of  preliminary  tests  prior  to  the  installation 
of  the  airfoil.  Flow  uniformity  was  verified  for  the  region  outside 
of  the  influence  of  the  four  test  section  boundary  layers  by  kiel  probe 
surveys.  The  longitudinal  static  pressure  gradient  [ (dPs/dx)/l/2pVco^] 
was  measured  to  be  0.01236  ft~l  (0.0406  m-1).  The  sidewall  boundary 
layer  at  the  balance  shaft  location  was  measured  at  several  Reynolds 
numbers;  the  resulting  data  are  presented  in  Figure  8. 
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The  force  measurements  phase  with  the  cantilevered  balance  was 
initiated  by  measuring  the  tare  forces  on  the  11.0  in  (279.4  mm)  diameter 
disk.  For  these  measurements  the  9.0  in  (228.6  mm)  airfoil  was  mounted 
on  the  opposite  wall  with  a  gap  of  0.004  (0.102  mm)  between  the  airfoil 
tip  and  the  disk.  The  results  are  listed  in  Table  I. 


TABLE  I 


Disk  Tare  Forces* 

Lift  and  Drag  Data  -  Cantilevered  Balance 
=  70.0  fps  (21.34  m/s) 

Re  =  330,000 

disk  diameter  =  11.0  in  (279.4  mm) 


a  (deg) 


c£  (±0.0014) 


9.34 

6.26 

3.21 

0.47 

-3.43 

-6.17 

-9.01 


0.0094 

0.0078 

0.0055 

0.0039 

0.0027 

-0.0008 

-0.0004 


cd  (±0.0004) 

0.0021 

0.0014 

0.0018 

0.0012 

0.0009 

0.0014 

0.0020 


*9.0  in  (228.6  mm)  airfoil  mounted  on  opposite  wall 
coefficients  nond linens ionali zed  by  c  =  9.0  in  (228.6  mm) 


It  should  be  noted  here  that  to  match  the  NACA  Reynolds  number 
conditions  (Re  =  330,000)  the  9.0  in  (228.6  mm)  airfoil  was  tested  at 
70  fps  (21.34  m/sec)  wheieas  the  18.0  in  (457.2  mm)  airfoil  was  tested 
at  37.0  fps  (11.28  m/sec). 

Lift  and  drag  were  next  measured  for  the  9.0  in  (228.6  mm)  airfoil 
while  operating  both  with  and  without  the  disk  and  for  the  18.0  in 
(457.2  mm)  airfoil.  The  resulting  c£  versus  a  data  (with  the  traditional 
corrections  applied)  which  were  judged  to  be  in  satisfactory  agreement 
with  the  NACA  reference  data  are  shown  in  Figure  3.  The  corresponding 
drag  polar  is  shown  in  Figure  4,  Here  the  measured  drag  data  are 
approximately  45.0  percent  high  at  upper  c£  values.  It  is  also  interesting 
to  observe  from  the  data  that  a  disk  on  the  balance-end  of  the  airfoil  is 
apparently  not  necessary. 
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At  this  time  the  proven  NACA  approach  of  using  separate  drag  measure 
ments  was  attempted.  The  drag  characteristics  of  the  9.0  in  (228.6  mm) 
airfoil  were  measured  by  the  downstream  wake  traverses  at  the  midspan 
of  the  airfoil.  The  drag  values  were  obtained  by  the  classical  procedure 
of  integrating  the  wake  momentum  deficits.  The  resulting  polar  from  the 
combination  of  the  balance-measured  lift  data  and  the  wake-measured  drag 
data  agreed  with  the  NACA  data,  Figure  5. 

What  was  the  source  of  error  in  the  balance-measured  drag  data? 

Three  possible  sources  were  considered:  (1)  end  gap  effects  between 
the  airfoil  tip  and  the  adjacent  sidewall,  (2)  drag  on  the  portion  of  the 
balance  shaft  between  the  sidewall  and  the  model,  and  (3)  contamination 
of  the  two-dimensional  flow  by  an  interaction  of  the  airfoil  and  the  side 
wall  boundary  layer.  The  effect  of  the  end  gap  between  the  airfoil  tip 
and  the  test  section  wall  was  first  investigated.  Lift  and  drag  data 
were  measured  at  a  constant  angle  of  attack  while  varying  the  end  gap 
from  0.001  in  to  0.010  in  (0.025  mm  to  0.254  mm).  For  this  range  of  end 
gap,  no  significant  change  in  c,  was  measured  and  less  than  a  1.0  percent 
change  in  c^  was  recorded.  Parkin  and  Kermeen  [11]  reported  similar 
results,  i.e.,  if  the  end  gap  is  sufficiently  small,  viscous  forces 
predominate  and  the  effect  of  the  end  gap  is  negligible. 

To  evaluate  the  drag  force  on  the  portion  of  the  balance  shaft 
exposed  to  the  flow  a  stub  spindle  was  fabricated.  The  stub  spindle 
was  mounted  in  the  cantilevered  balance  and  extended  0.002  in  (0.508  mm) 
into  the  flow.  This  was  the  typical  operating  clearance  at  the  balance 
end  of  the  airfoil.  The  9.0  in  (228.6  mm)  airfoil  was  mounted  on  the 
opposite  wall  and  was  maintained  at  a  minimum  distance  from  the  spindle. 
Under  these  conditions,  the  flow  in  the  vicinity  of  the  model-sidewall 
intersection  is  closely  duplicated  and  the  balance  measured  only  the 
forces  on  the  stub  spindle.  The  resulting  data  are  presented  in 
Table  II. 


TABLE  II 

Stub  Spindle  Drag  Measurements* 
=  70.0  fps  (21.34  m/sec) 

Re*  =  330,000 


a  (deg) 

c,  (+0.0004) 
d 

0.00 

0.0000 

7.10 

0.0002 

7.77 

0.0006 

8.90 

0.0008 

10.80 

0.0013 

*9.0  in  (228.6  mm)  airfoil  mounted  on  opposite  wall 
coefficient  nondimensionalized  by  c  =  9.0  in  (228.6  mm) 
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Clearly  the  supporting  shaft  was  not  responsible  for  the  high  measured 
values  of  c^. 

Thus,  only  the  contamination  of  the  two-dimensional  flow  by  the 
interaction  of  the  airfoil  with  the  sidewall  boundary  layer  remained 
as  a  postulated  cause  of  the  erroneous  drag  measurement.  Evidence  of 
such  an  interaction  was  observed  during  the  wake  measurement  phase  where 
a  secondary  wake  was  measured  near  the  sidewall.  Figure  9.  As  previously 
discussed  the  removal  of  the  sidewall  boundary  layer  is  not  practical  in 
water  tunnel  applications  so  that  the  alternate  approach  of  developing  a 
correction  procedure  was  chosen. 

DEVELOPMENT  OF  THE  CORRECTION  PROCEDURE 

As  shown  in  Figure  1,  2  and  4  a  similar  variation  between  balance- 
measured  drag  data  and  the  NACA  two-dimensional  section  characteristics 
existed  in  Kermeen’s,  Daily's  and  the  ARL/PSU  data.  Because  this 
difference  is  typical  only  of  balance-measured  data  in  which  the  entire 
force  exerted  or  the  airfoil/hydrofoil  is  measured,  it  was  assumed  that 
this  increment  in  drag,  Ac  ,  was  due  to  three-dimensional  flow  effects 
on  the  model.  Such  three-dimensional  flow  effects  can  be  generated  when 
a  strut  intersects  a  flat  surface  in  the  presence  of  a  nonuniform  flow. 

The  resulting  secondary  flow  -  the  so-called  horseshoe  vortex,  Figure  10  - 
engulfs  the  strut-wall  intersection  and  produces  a  region  of  contaminated 
two-dimensional  flow.  This  type  of  secondary  flow  can  be  generated  in 
airfoil/hydrofoil  testing  when  an  airfoil  or  hydrofoil  that  spans  the 
test  section  intersects  the  test  section  wall  in  the  presence  of  the 
sidewall  boundary  layer. 

A  study  of  this  problem  using  flow  visualization  techniques  was 
recently  completed  by  Barber  [6]  in  which  he  investigated  the  additional 
drag  that  is  created  by  a  strut  protruding  from  a  wall  as  a  function  of 
the  incoming  boundary  layer  thicknesses.  He  found  that  the  size  of  the 
horseshoe  vortex  varied  directly  with  the  thickenss  of  the  incoming 
boundary  layer.  He  also  found  that  the  portion  of  the  airfoil  where 
flow  separation  occurred  varied  inversely  with  the  size  of  the  horse¬ 
shoe  vortex.  He  concluded  that  with  a  large  horseshoe  vortex,  viscous 
effects  caused  high  energy  fluid  to  be  entrained  in  the  corner  where  the 
airfoil  trailing  edge  and  wall  intersects  as  shown  in  Figure  11.  This 
influx  of  high  energy  fluid  enables  the  flow  to  withstand  better  the 
adverse  pressure  gradient  existing  in  the  corner  and,  consequently, 
retards  flow  separation.  As  illustrated  in  Figure  11,  a  thin  vortex 
is  not  able  to  entrain  as  much  of  the  high  energy  fluid  and  a  larger 
separated  zone  exists. 

Hawthorne  [11]  derived  the  following  expression  for  the  energy  in 
secondary  flows,  Dg*  created  by  strut-wall  intersections: 


D 

e 


144  V  2  c2  (t/c)4  f (n) 

GO 

25  [1  +  (l/2)(t/c)2] 


(1) 


?*«*  ir-  •****.^9*  iv 


* 
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where 


f(n)  = 


2  , 
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2  TT 
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2_  ilLzl 

4n  L2+v 


2  i  2 

+  -  ~n2^  log  n  + 
(1+n)  e 


1+n 


4n 


and 

n  =  4  [1  +  (1/2)  (t/ c)  ]/(15ir  <5*/c) 


(2) 


(3) 


Hawthorne's  relationship  between  f(n)  and  boundary  layer  displacement 
thickness  (6*/c)  is  shown  in  Figure  12.  These  data  are  for  a  bicusped 
strut  profile  in  an  exponential  boundary  layer  with  strut  thickenss-to- chord 
ratios  of  .05  and  .25.  Hawthorne's  figure  shows  that  f(n)  increases  with 
5 */c  to  a  maximum  value  at  <5*/c  *  0.1.  Equation  (1)  states  that  the 
energy  in  these  secondary  flows  is  proportional  to  airfoil  thickness  to 
the  fourth  power  and  reaches  a  maximum  when  6  /c  is  approximately  0.1. 
Although  the  theory  does  not  hold  for  all  airfoil  shapes  or  boundary 
layer  profiles,  it  is  probably  fair  to  assume  that  generally  the  energy 
in  secondary  flows  for  this  type  of  airfoil-tunnel  wall  intersection 
is 


* 

D 

e 


(t/c)4  f (n) 


(4) 


If  the  function  fj^n)  in  Figure  12  is  linearized  over  the  portion  of 
the  curve  0.0  <  6  /c  <  0.1,  then 


f  (n) 


f(n)  /(0.D 

IT)  3.X 


(5  /c)  =  K2  (6/c) 


(5) 


In  Equation  (5),  the  momentum  thickness,  6  ,  has  been  assumed  proportional 
to  the  more  frequently  documented  boundary  layer  thickness,  6.  Thus,  Dg* 
becomes 


De*  =  K3  (t/c)4  (6/c) 


(6) 


Functionally,  the  drag  correction  was  assumed  to  take  the  following 

form: 


Acd 


81 


'V 


cd,  a,  AR) 


(7) 


The  inverse  relationship  between  Ac.  and  D  has  been  established  by 
Barber  [6].  Therefore,  Equation  (7)  can  be  written  as 
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g2  (c  ,  cd,  a,  AR) 

Ac  ,  =  K.  - - - 

(6/c)(t/c) 


(8) 


The  effects  of  c^,  c^,  a  and  AR  were  derived  empirically  from  the  available 
experimental  data. 


The  required  Ac^  corrections  for  the  three  sets  of  drag  data  are 
shown  by  the  dash-dot  curves  in  Figures  1,  2  and  4.  The  three  sets  of 
data  have  several  similarities.  The  deviation  of  drag  is  zero  at  c^=0.0 
for  all  three  investigations.  The  Ac^  curves  for  all  three  experiments 
increase  to  the  points  where  c^  is  no  longer  constant  and  then  decrease. 

If  c^  represents  the  slope  ofathe  lift  curve  in  the  linear  portion  of 
a 

the  c^  versus  a  curve,  then  the  shape  of  the  Ac^  versus  a  curve  seems  to 

vary  as  [co  /c^  ]^/2.  for  all  three  experiments.  It  was  assumed 

a  aQ 

that  this  slope  variation  represented  the  angle  of  attack,  a,  contribution 
to  the  drag  correction  so  that  Ac^  can  be  expressed  as 


Ac,  =  Kc 
a  o 


a  a. 


1/2 


g3(V  Cd’  AR) 


(5/ c) ( t / c) 


(9) 


As  the  experimental  data  from  the  three  studies  show,  Ac^  increases 
directly  with  c^;  and  Ac^  is  zero  at  c^  =  0.0.  This  also  implies  that 
the  balances  measure  the  correct  value  of  c^  at  c^  =  0.0,  namely,  c^  . 

Thus,  the  c,  was  included  as  the  c,  term  which  "individualizes"  the° 
d  a 

correction  procedure  to  specific  airfoils.  When  the  linear  dependence 
on  c^  and  the  c^  term  are  introduced  in  Equation  (9),  Ac^  becomes 


cd. 


1/2 


(6/c) (t/c)^ 


g4(AR) 


(10) 


The  effect  of  aspect  ratio  was  approximated  by 

g4(AR)  =  (AR)b  .  (11) 

The  values  of  K.  and  b  were  determined  empirically  from  the  experimental 
data.  It  was  found  that  the  best  fit  to  the  three  sets  of  data  was  obtained 
for  K,  =  1.5  x  10~^  and  b  =  -1/2.  With  the  evaluation  of  these  two  constants 
the  final  form  of  the  equation  to  correct  the  balance  measured  drag  data  for 
the  effect  of  the  sidewall  boundary  layer  is 
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Acd 


=  1.5  x  10 


(cit)(cdo)(cVC*a  )1/2 
(6/c)(t/c)4(AR)1/2 


(12) 


In  summary,  the  proposed  correction  to  balance-measured  drag  data 

for  the  effect  of  the  sidewall  boundary  layer  is  a  function  of  the  airfoil/ 

hydrofoil  geometry  (t/c  and  AR),  the  thickness  of  the  sidewall  boundary 

layer  (6/c),  and  the  accurately  measured  balance  data  (c,  and  c„  versus  a). 

do  Ji 

APPLICATION  AND  DISCUSSION 


Shown  in  Figures  13,  14  and  15  are  the  results  of  applying  the  correction 
procedure  to  the  Kermeen,  Daily  and  ARL/PSU  data,  respectively.  For  Daily's 
half-span  measurements  when  the  balance  sensed  the  effect  of  only  one 
sidewall  boundary  layer,  the  total  AR  value  was  used  but  only  one  half  of 
the  computed  Ac^  value  from  Equation  (12)  was  applied  to  the  measured  data. 

As  can  be  seen  in  these  three  figures  the  corrected  drag  polars  are  in 
relatively  good  agreement  with  the  accepted  reference  data. 

In  practice,  this  correction  for  the  sidewall  boundary  layer. 

Equation  (12),  is  applied  to  the  drag  data  after  the  tare  readings  and 
the  traditional  corrections  have  been  applied.  The  technique  will  be 
illustrated  by  applying  the  procedure  to  data  measured  by  Ward  [8]  at 
CIT  for  the  Canadian  Defense  Research  Establishment  Atlantic  (DREA) . 

These  data  were  located  by  the  authors  after  the  completion  of  Jacobs' 
initial  studies  and  were  not  included  in  the  development  of  Equation  (12). 

Ward  used  a  three  component,  cantilevered,  mechanical  balance  to 
measure  lift,  drag  and  pitching  moment  on  a  6.0  in  by  6.0  in  (152.4  mm 
by  152.4  mm)  NACA  16.309  hydrofoil  in  the  CIT  High  Speed  Water  Tunnel. 

The  resulting  noncavitating  data  at  50.0  fps  (15.24  m/s)  with  the  tare 
corrections  included  are  shown  by  the  open  circles  in  Figures  16  and  17. 

Also  shown  as  solid  lines  in  these  figures  are  the  NACA  reference  data 
(after  correction  for  compressability  effects)  measured  by  Lindsey,  et.al. 

(12]  at  a  Mach  number  of  0.3.  The  data  were  further  corrected  according 
to  Pope  [13]  for  solid  blockage,  wake  blockage  and  lift  effect  (stream¬ 
line  curvature).  The  sidewall  of  the  test  section  was  adjusted  to  eliminate 
the  horizontal  buoyancy  effect.  The  results  of  applying  these  traditional 
corrections  to  the  data  are  shown  by  the  open  squares  in  Figures  16  and  17. 

The  required  Ac^  correction  for  Wards'  data  is  shown  by  the  dash-dot 
curve  in  Figure  17.  Again  it  is  interesting  to  note  that  the  balance 
has  measured  the  correct  drag  value  at  the  zero  lift  condition. 

From  Reference  [8]  the  following  parameters  were  obtained  for  the 
application  of  Equation  (12): 

cd  =»  0.0009  AR  =  1.00  (t/c)  =  0.09 

o 
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The  term  is,  of  course,  the  corrected  lift  coefficient  (open  squares) 
at  each  a.  The  (c^  /c^  51/2  term  was  computed  by  fitting  a  differentiable 

a  a 

o 

mathematical  spline  curve  through  the  corrected  c^  versus  a  data.  Ward  [14] 
documents  the  test  section  boundary  layer  characteristics  at  the  balance 
shaft  location  of  the  HSWT.  The  average  value  of  6/c  at  50.0  fps  (15.24  m/s) 
is  0.108.  With  these  values  Equation  (12)  can  now  be  evaluated.  The  resulting 
Ac^  values  are  shown  as  the  solid  squares  in  Figure  17.  When  these  computed 
Ac^  values  are  applied  to  Ward's  drag  polar  the  corrected  data  are  shown  by 
the  solid  circles.  The  agreement  is  certainly  encouraging. 

CONCLUSIONS,  LIMITATIONS  AND  RECOMMENDATIONS 

The  results  of  this  ARL/PSU  experimental  investigation  and  subsequent 
data  analysis  have  revealed  or  reaffirmed  several  important  conclusions 
relative  to  balance-oriented,  two-dimensional  airf oil/hydrofoil  testing: 

(1)  The  effect  of  a  small  end  gap  between  the  airfoil  tip  and  the  channel 
wall  is  negligible  provided  the  gap-to-chord  ratio  is  _<  0.002. 

(2)  It  was  also  demonstrated  that  the  use  of  a  flush-mounted  disk  at  the 
balance-end  of  the  airfoil  is  not  required,  at  least  in  noncavitating  flows. 

(3)  In  the  absence  of  the  disk  and  for  gap-to-chord  ratios  £  0.001  the 
effects  of  flow  in  the  region  of  the  supporting  shaft  are  negligible. 

(4)  With  the  application  of  only  traditional  and  tare  corrections,  drag 
polars  in  agreement  with  accepted  reference  data  can  be  obtained 

by  combining  the  balance-measured  c^  values  and  c^  data  from  wake 
traverses. 

(5)  The  disagreement  between  balance-measured  drag  data  and  the  accepted 
reference  values  is  primarily  the  result  of  the  interaction  of  the 
airfoil  and  the  sidewall  boundary  layer. 

(6)  The  effect  of  the  sidewall  boundary  layer  on  balance  measured  drag 
data  can  be  accounted  for  by  the  application  of  Equation  (12) . 

The  previous  conclusions  are  not  without  some  limitations.  The 
empirical  development  of  Equation  (12)  was  conducted  in  the  absence  of 
data  obtained  from  studies  in  which  there  was  a  significant  variation 
in  aspect  ratio  or  thickness  to  chord  ratio.  The  linearized  adaptation 
of  Hawthornes  f(n)  -  curve  is  only  valid  for  6*/c  £  0.1.  For  values  of 
<5*/c  >0.1  a  different  approximation  of  f(n)  would  be  required.  Equation 
(12)  has  been  developed  in  terms  of  the  boundary  layer  thickness  measured 
at  the  balance  shaft  location.  Since  all  of  the  subject  foils  were  mounted 
at  their  midchord  position  it  is  not  totally  apparent  that  the  midchord 
location  is  the  appropriate  position  for  the  boundary  layer  measurements. 


Obviously  one  of  the  recommendations  for  future  study  would  be 
the  acquisition  of  a  larger  data  base,  particularly  with  respect  to  AR 
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and  t/c  variation,  for  further  refinement  of  Equation  (12).  The  effects 
of  the  sidewall  boundary  layer  should  be  further  investigated  by  wind 
tunnel  tests  in  which  the  horseshoe  vortex  can  be  removed  by  suction  or 
blowing.  And,  of  course,  an  attempt  should  be  made  to  extend  the  proposed 
correction  procedure  to  hydrofoils  operating  in  the  cavitating  flow 
regimes. 
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Figuru  1.  Kermeen's  drag  polar  compared  with  NACA  published  data  for  a  NACA  4412 
airfoil.  Also  shown  is  the  value  of  AC^  required  to  correct  Kermeen's 
curve  to  the  NACA  curve. 
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Figure-  2.  Daily’s  drag  polar  compared  with  NACA  published  data  for  a  NACA  4412 
airfoil.  Also  shown  is  the  value  of  AC^  required  to  correct  Daily's 
curve  to  the  NACA  curve. 


V 


ARL/PSU  drag  polar  (traditional  corrections  only  applied)  for  c  =  9.0 
(228.6  mm)  and  c  =  18.0  (457.2  mm)  airfoils  without  disks.  Also  shown 
the  value  of  AC,  required  to  correct  the  ARL/PSU  data  to  the  NACA  curv- 


Comparison  of  NACA  data  with  ARL/PSU  and  data  measured  by 
different  techniques. 
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APPLIEO  RESEARCH  LABORATORY,  FLUIDS  ENGINEERING  DEPARTMENT 
THE  PENNSYLVANIA  STATE  UNIVERSITY 

P.  0.  BOX  30,  STATE  COLLEGE,  PA  16801 _ (814)  865-1741 


SUBSONIC  WIND  TUNNEL 


(USA) 


DESCRIPTION  0?  FACILITY: 


Closed  Circuit 


TYPE  OF  DRIVE  SYSTEM: 


Axial~Flow  Blower  Variable-Speed 


TOTAL  MOTOR  POWER: 


.150  HP  (111.8  kv) 


WORKING  SECTION  MAX.  VELOCITY:  45.72  a/s 


INSTRUMENTATION:  Automatic  scanning  and  rotating  mechanisms,  pressure 
sensors,  hot-wires 


TURBULENCE  LEVEL:  0.2  percent 


PROPELLER  OR  MODEL  SIZE  RANGE:  Up  to  635.0  mm  dia. 


TESTS  PERFORMED:  Pressure  distributions  over  hydrodynamic  shapes.  Velocity 
profiles  in  propeller  planes.  Hot-wire  measurements.  Wall  interaction,  effects 
on  hydrofoils. 


OTHER  REMARKS:  The  tunnel  is  subsonic.  It  is  used  for  basic  and  applied 
research.  It  is  well  instrumented  for  the  measurement  of  model  boundary 
layers  and  their  wake  turbulence. 


PUBLISHED  DESCRIPTION:  ARL/PSU  Report  NORD  16597-56,  Lehman,  1959 


Figure  6.  ARL/PSU  subsonic  wind  tunnel. 
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Figure  9.  Wake  profiles  in  the  vicinity  of  the  airfoil-disk  intersection 
shoving  the  growth  of  a  second  wake  behind  airfoil  upper 
surface. 
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Figure  15-  ARL/PSU  results  corrected  by  Equation  (12). 
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